Unified theory of phase separation and charge ordering in doped manganite 
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A unified theory is developed to explain various types of electronic collective behaviors in doped 
manganites Ri-^jX^iMnOs (R = La, Pr,Nd etc. and X = Ca, Sr, Ba etc.). Starting from a realistic 
electronic model, we derive an effective Hamiltonianis by ultilizing the projection perturbation 
techniques and develop a spin-charge-orbital coherent state theory, in which the Jahn- Teller effect 
and the orbital degeneracy of eg electrons in Mn ions are taken into account. Physically, the 
experimentally observed charge ordering state and electronic phase separation are two macroscopic 
quantum phenomena with opposite physical mechanisms, and their physical origins are elucidated in 
this theory. Interplay of the Jahn- Teller effect, the lattice distortion as well as the double exchange 
mechanism leads to different magnetic structures and to different charge ordering patterns and phase 
separation. 

PACS numbers; 75.30.Mb, 71.27.-|-a, 75.10.-b, and 75.45-hj 



I. INTRODUCTION 

Charge jQpdering states (CO)ErLl and electronic phase 
separationBiJ (PS) are two of the macroscopic quantum 
phenomena observed experimentally in doped mangan- 
ites Ri_a;X:rMn03 (R=La, Pr, Nd etc. and X=Ca, Sr, 
Ba etc.). This seems puzzling since these two phenom- 
ena have completely opposite physical mechanisms. CO 
near a; = 0.5 is a regular alignment of Mn'^+ and Mn^+ in 
the real space. It is well known that the Wigner lattice is 
expected to be stabilized when the repulsive potential be- 
tween charge carriers dominates over the kinetic energy 
of the carriers. In this respect CO is expected to form in 
manganites due to strong repulsion between charge car- 
riers. Oppositely, PS in doped manganites near x — 
is characteristic of two regions of rich- and poor-density 
of charge carriers with ferromagnetic (FM) and antifer- 
romagnetic (AF) correlations, respectively. A uniform- 
density state is unstable when the charge carriers are 
subjected to a strong attractive interaction, as discussed 
in high Tc superconductors.cHl It should be an strong at- 
traction which drives the charge carriers to the electronic 
PS. On the other hand, various types of magnetic struc- 
tures and orbital ordering states were also observed ex- 
perimentally. CO and PS are definitely associated with 
these structures. For instance, CO with the (tt, tf, 0) pat- 
terns occurs under the C-type antiferromagnetic (AF) 
background, and PS near a; = occurs under the A-type 
AF background. The field-induced melting effect of CO 
shows that the CO decreases and eventually disappears 
while an external magnetic field increases.EJ The field de- 
stroys the AF correlation, and the disappearance of CO 
and AF indicates their close relations and the possible re- 
lation between AF correlation and repulsive interaction 
of charge carriers. Hence experimental observations of 
CO and PS in manganites strongly suggest that the sign 
of effective interaction between the charge carriers should 
depend on the dopant concentrations and the magnetic 
structures. 



There has been considerable theoretical work moti- 
vated by the experimental research on manganese oxides. 
Most theoretical efforts focus on understanding metallic 
ferromagnetiain and its connection to unusual transport 



properties 



The scenario of double exchange mecha- 



nism is extensively accepted to explain the metallic ferro- 
magnetism. However, we still lack a comprehensive pic- 
ture for the physical origins of the PS and CO in doped 
manganites. In a simplified one-bajiKjl -model, PS was 
studied analytically and numerically.Eao An attraction 
between the charge carriers caused by the superexchange 
coupling is rosponsible for the instability of a uniform- 
density state.ll3 In the vicinity of a; = 0.5, the mechanisms 
of both long-range Coulomb interaction and|-tl|ie,particle- 
hole interaction for the CO were proposed.EJEiicS How- 
ever, PS and CO cannot be explained simultaneously in 
the same one-band model as strong long-range Coulomb 
interaction does not favor forming the PS near x = 0. 
Other properties, such as the ferromagnetism at x = 
and anomaly optical conductivity, also suggest that the 
double degeneracy of eg orbital, which is neglected in 
the one-band model, should be included. The phase 
diagrams of doped manganites, especially for r«iagnetic 
and orbital ordering, have been investigated.E3c3_Recent 
overviews for doped manganites are seen in Ref.E3c3. 

In this paper, we explore the origins of PS and CO 
in doped manganites and establish a unified theory for 
these two phenomena. The paper is organized as follows. 
An effective Hamiltonian is derived in Section II. Start- 
ing from a realistic electronic Hamiltonian with strong 
electron correlations, several virtual processes of superex- 
change are considered and an effective Hamiltonian is de- 
rived by means of the projective perturbation approach 
and Schwinger boson formalism. A theory of spin-charge- 
orbital coherent state is presented in Section III. Close 
connections of the PS and the CO with various types of 
AF are elucidated in Section IV. We also show that the 
Jahn- Teller (JT) effect and lattice distortion play impor- 
tant roles in stabilizing the magnetic structures. Some 



discussions and a brief summary are given in Section V. 
Detailed derivation of the effective Hamiltonian up to the 
second order is presented in Appendix A, and the three- 
site interaction terms are given in Appendix B. 



II. EFFECTIVE HAMILTONIAN: A 
PROJECTION PERTURBATION APPROACH 

A. Model Hamiltonian 

Doped manganese oxides, Ri_a;Xa;Mn03, can be re- 
garded as a mixture of Mn'^'*' (3d"') and Mn"'+ (3d"^) ions. 
The three electrons in the outer shell of Mn'*+ are almost 
localized in i2g orbit to form spin maximal state with 
S = 3/2. In the manganese ion, Mn'^+, apart from the 
three localizes electrons in t2g orbit, the fourth d electron 
locates at eg orbit which is doubly degenerated. Cg elec- 
trons can become delocalized with increasing x. From 
the above reasoning, an electronic Hamiltonian with or- 
bital degeneracy is put forward to describe the dominant 
electron-electron interaction of the system 
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are the creation and annihila- 
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tion operators of Cg electron at the orbit 7 
where and \z) oc {Sz"^ — r^)/v3 and \z) oc x^ ~ y^ , re- 
spectively.) of site i with spin a {—1, J,), respectively. 
Si, 7 is the spin operator of eg elec- 
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Si is the total maximal spin operator of the 
three t2g electrons. Here the transfer integrals in the 
model are assumed to take a Slater-Koster form given by 
the hybridization between e^ orbit and nearest oxygen p 
orbit, and the model has been extensively studiec l to u p-, 
derstand physics of doped manganite provskites.E^lEa^EII 
It was already realized that the Hund's rule coupling, 3h, 
between e^ electron and three t2g electrons is very strong, 
which is the main originjOf metallic ferromagnetism in the 
range of 0.1 < x < 0.5£3 Most eg electrons favor form- 
ing S — 2 spins with localized t2g electrons in Mn^+. 
Usually only the strong Hund's rule coupling 3h is taken 
into account. However, the on-site Coulomb interaction 
U (intra orbit 7 = 7') and U' (inter orbits 7 ^ 7') are 
also dominant energy scales and usually larger than J^S. 
We will see that the strong on-site correlation play an es- 
sential role of the electronic collective behaviors. 
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FIG. 1. Five processes that generate the effective Hamilto- 
nian along the z-axis. (a) a direct hoping of Sg electron from 
one site to its nearest neighbour site. The spin of electron 
must be parallel to the localised spin at the same site. The 
(b-e) are four virtual ones to the superexchange couplings 

(b) the mediate state is a single occupancy with spin 5'— 1/2 

(c) the mediate state is a double occupancy on the same orbit 
(d).the mediate state is a double occupancy with spin 5-1-1 
on different orbits; and (e) the mediate state is a double oc- 
cupancy with spin S on different orbits. These five processes 
corresponds to the five terms in Eq.(2), respectively. The in- 
teractions along other x- and y-axises are derived by using 
the symmetry of rotation. 

Besides the part Hamiltonian of the conduction elec- 
trons, we need take into account other parts of interac- 
tion which are believed to affect the phase diagram of 
the doped manganites. First, a tiny hopping between t2g 
electrons produces a weak superexchange AF coupling 
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Apart from the electronic part of interactions, the JX 
effect leads to a static distortion and mixes eg orbits.E^ 
According to KanamoriJ^J at x=0, the primary lattice 
distortion is a staggered (tt, tt, tt) tetragonal distortion of 
the oxygen octahedra surrounding the Mn sites, driven 
by a Jahn- Teller splitting of the outer Mn d levels. 
Kanamori's deduction was subsequently confirmed by 



more detailed studies of the structureHa Hence an ef- 
fective Hamiltonian is introduced phenomenologically for 
the JT effect, 

HjT = k Y, (Ni,a,^-Ni,a,<T)(Nj-„,^--Nj,a^^O- 

where a depends on the direction of r,; — Vj . The tetrag- 
onal crystal field H^ wilLLcad to the anisotropic magnetic 
structure of the systcm:E3 
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This two parts of interaction is independent of spin. 

In short, combining the JT effect and the lattice distor- 
tion caused by the crystal field, the total effective Hamil- 
tonian is found to be 
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B. A projection perturbation approach 

In the large Jh limit and when the number of electrons 
is not greater than the number of lattice sites, each site 
is occupied by at most one Cg electron. And the electron 
must form a spin S +1/2 state with the localized spin on 
the same site. The process can be realized by introducing 
the projection operator 



p-Uip^'^ + T^p^s)- 
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For any state \a), P\a) is the component of the state 
with holes and single occupancies with spin S+1/2. This 
technique has been applied eytensively to study the one- 
band Kondo lattice modelEj^Ell The projection operators 
are defined by 
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Hence, in the infinite Jh limit, the Hamiltonian in Eq. 
M is reduced to 



Hi"^ = PH,P 



(4) 



as shown in Fig. 1(a). This is the double exchange 
model with orbital degeneracy|_-. With the help of the 
Schrieffer- Wolff transformation,L3 the finite, but large Jh 
effect is taken into account by superexchange processes in 
the second-order projective perturbation approach: (1) 



Mn3+Mn4+ ^ Mn^+Mn^+as shown in Fig. 1(b) and (2) 
Mn3+Mn3+ <^ Mn'*+Mn2+ as shown in Figs.l(c - e). Af- 
ter considering the second order perturbation correction 
to Eq. (0), we obtain 
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where AEa are the energy difference of the mediate state 
and the initial state as shown in Figs, l(b-e). Qa 
(a — &, c, d, e) are the projection operators for the in- 
termediate states. The derivation of Eq. (0) and the 
explicit expressions of A_Eq and Qa are presented in Ap- 
pendix. 



C. Effective Hamiltonian in Schwinger boson 
representation. 

To simplify the notations, we expresS|-the Hamiltonian 
in the Schwinger-boson representation.LlI The represen- 
tation was introduced, to describe the one-band double 
exchange by Sarker,L3 and here we generalize it to the 
our model by introducing another type of boson for or- 
bital degree of freedom. Define 



p+ c^. p+ = L 






where a] ^ and Ui^cr are the Schwinger boson operators 
for spin: 



I'^j,],) 






b] Q, and bi^a are the Schwinger boson operators for orbital 



degrees of freedom with a 
direction of r^j. b} ^\0) = 



X, y, z, which depends on the 
\z) and big |0) = \z) . Other 

two components are not independent and related to the 

z-component by a transformation: 
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fj and fi are fermion operators for charge carrier: rii — 
fl fi = 1 means that there is one charge e on the site, 
i.e., Mn'^+, and m = means Mn'*+. The constraints for 
the Schwinger boson and fermions are 



'^It'^^'T + '^li^'^i = 2S + \ 



In the representation, the effective Hamiltonian for con- 
duction electrons is written as, 



He - Hi"^ + Hi'''^ + Hi^^ + Hi"-^ + Hi" 



(6) 



Each term corresponds to one of the processes shown in 
Fig. 1, and is expressed as 

iy('')-_V ^ flt a b^ b- fU 
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where Si is a spin operator with S", and Si is a spin oper- 
ator with S -f 1/2 as a FM combination of the locahzed 
spin and itinerant electron at the same site. The opera- 
tors P are the projection operators for charge and orbits: 

P'^ = n^n" 

Plh = l~ Ui, 

P — n" 

where n" = b\jji^a- Finally, it is worth mentioning 
that we just keep the two-site interactions and neglect 
three-site interactions in Eq. (0). The three-site terms 
describe indirect hopping processes between the next- 
nearest-neighbor site via the intermediate states and is- 
believed to be relatively small compared with the direct 
hopping terms. 

The other terms in Eq. (g) become 

Haf = Jaf 2_^ \^i ' ^j ^ ^ ) ' 

HjT^kY,K-nf){n'^-nf), 

i 

respectively, where Si = Si(l — Jii) + ^'^.^ Sirii. 

Hence up to the order of I^/Jh the total effec- 
tive Hamiltonian in the representation of the Schwinger 
bosons for both spin and orbit is 



Heff - He 
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Here each term should be restricted in the projected 
space. Approximately, -ffe// in Eq- (Q) ^^id Htotai in 
Eq. (||) are expected to describe the same low energy 
physics in the large 3h case. 

The present theory is based on Eq. (R), in terms of 
which we are able to establish a unified description for 
the electronic behaviors in doped manganites for the first 
time. In this paper, the model parameters are roughly 
estimated from the excitation euergies of Mn ions and the 
functional density calculations:cj we take t = OAleV as 
energy unit. 2t/JH{2S+l) = 0.35; t/{U + JHS) ^ 0.042; 
t/{U' + I J + JhS) = 0.056; t/{U' - i J) = 0.106; and 
Jaf — 0.001. All the phase diagrams in this paper are 
established on this set of parameters. 



III. SPIN-CHARGE-ORBITAL COHERENT 
STATE FORMALISM 

The second order projection perturbation approach in- 
cludes part of strong correlation between conduction elec- 
trons and localized spin, and removes the direct Hund 
couphng 3h- Some properties of the model [Eq. (Q)] 
become clearer than the original one [Eq. (|l|)]. For in- 
stance. He describes the double exchange mechanism 
for ferromagnetism; He describes a particle-hole inter- 
action with an AF coupling; He and He describe the 
AF superexchange couplings with the same and differ- 
ent orbitals, respectively, and He -^ describes a FM su- 
perexchange coupling with different orbitals. Each term 
becomes predominate at some point in the phase space. 
However, the effective Hamiltonian [Eq. (^] still seems 
to be very complicated, and it is still very hard to fix 
its physical properties. From the study of the famous 
t-J model, which is derived from the one-band Hubbard 
model by the same projection perturbation, we learn that 
this is just a first step to understand physics of the elec- 
tronic Hamiltonian in Eq. (yj) in the strong Hund cou- 
pling case. 



A. Spin-charge-orbital coherent state 

To investigate the effective Hamiltonian, we apply the 
spin-coherent state mean-field theory. We introduce two 
polar parameters 9i and (j)i for spin bosons at the site i, 
and two pararpeters a^ and Pi for orbital bosons. Follow- 
ing Auerbach,c3 we define the spin-charge-orbital coher- 
ent state 

\0i,(l>i,ai,/3i,^i) = \Oi,(j>i)g\0) 

+ |e.,</'»)s+i/2«)|a^,A)i/2e^//|0) 

where 
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and ^ is an anticonimuting Grassmann variable. Define 
The Hamiltonian function is 

($1$) 
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ni^a = COS ( -;r + Oa 1 — sin 2(3q sin a^ sin — , 

"■I, a = sin^ \~^ ^ ^<^) ~'~ ^^"^ ^'^" ^^^ '^^ ^^^^ "T 
with 

cos Qij — cos 61^ cos 9j + sin ffj sin 0j cos{(pi — (/)j) 
and Sx — —T^/3, Sy — 7r/3, and Sz = 0. 

B. Mean Field Approximation 

Ti. includes the fourth powers of Grassmann variables. 
These terms are hard to integrate, and an approximation 
is needed. 

///./]/, « (///.) /]/, + ///. (/j/,) - (fin) (/J/,) • 



fjfi) is taken to the density of electrons in a density- 
uniform state, and of sublattice-depcndence when we 
consider a charge ordering state. The polar parame- 
ters are alsp-|treated in the mean field theory, following 
de Gennes.cSl The angles between two nearest-neighbor 
spins are taken to be (— l)*''"*~""^)'^(5cf If all Qa are equal 
to zero, it is a FM state. For a A-type AF, Qa are zero 
among the x-y plane, and Qz is non-zero along the z-axis. 
It is worth stressing that in this paper the antifcrromag- 
netism along a specific direction does not mean that Qa 
must be tt. For a C-type AF, Qz is taken to be zero 
and Qx,y are non-zero. G-type AF means that all Qa 
are equal, but non-zero. The parameters for orbital de- 
grees of freedom are treated similarly as those for spins. 
In this paper, we limit our discussion only in the case of 
T = 0, i.e., the ground state. The phase diagrams are 
established by minimizing the free energy. 



IV. CHARGE ORDERING AND MAGNETIC 
STRUCTURES 

Using the spin-charge-orbital coherent state theory, we 
focus on three phenomena which occur at different den- 
sities of doping: (a) ferromagnetism and A-type antifer- 
romagnetism at x = 0; (b) phase separation at small 
doping; and (c) charge ordering aX x = 0.5. 



A. Ferromagnetisni at s = 

The magnetic structure of the parent compound {x = 
0) has also very decisive impact on the electronic prop- 
erties near the point, and is sensitive to the model pa- 
rameters. The effective spin superexchange coupling is 
approximately 



Jeff/t 



U + JhS[S+1/2Y 
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which depends not only on the model parameters, but 
also the orbital orderings. It is worth noticing that the 
factor before {Pfj) is always negative, and the factor be- 
fore {Pij) is positive. When U is taken to be infinite, 
the factor before (-P/,) vanishes. In other words, the AF 
coupling in He is suppressed completely. As the FM 
coupling in He is always stronger than the AF coupling 
in Hi (i.e., the factor before {P^j) is always negative), 
the FM coupling becomes predominant, and the ground 
state becomes FM at low temperatures if we do not take 
into account the tiny AF coupling of J^i?0 In the case, 
the magnetic structure is independent of the JT effect 
and the lattice distortion, i.e., the orbital distribution of 
conduction electrons. When U is finite, the strength k of 
the JT effect and the lattice distortion coefficient e^ affect 
the magnetic structure by adjusting orbital ordering: the 
JT effect favors to form a "Neel's type AF" orbital-orbital 
correlation, which enhances the FM coupling through the 
process (e) in Fig. 1(e), while the lattice distortion tends 
to force the orbital to polarize, which increases the AF 
coupling through the processes ( b - d ) in Fig. l(b-d). 
The phase diagram in Fig. 2 shows that the FM survives 
at finite U for a small e^ and large k and evolves into A- 
type AF for a large ez and small k (We choose the model 
parameters in Section III.E.) Hence the ferromagnetism 
at a: = comes from the superexchange process in Fig. 1 
(e), not the double exchange mechanism in Fig. 1 (a) as 
the direct hopping is prohibited due to the strong cou- 
pling at a; = 0. The anisotropy of magnetism at x=0 orig- 
inates from the crystal field along the z-direction, which 
forces orbital degrees to form a "FM"-like state along 
that direction while the "AF" ordering remains among 
the z-y planes. It is worth mentioning that the ferromag- 
netic superexchange mechanism appears only when the 
orbital degeneracy of Cg electrons is taken into account. 
For a simplified one-band ferromagnetic Kondo lattice 
model, the ground state at half filling is always antiferro- 
magnetic for any finite Zh and UJ13 Hence to understand 
the ferromagnetism, we have to take the orbital degener- 
acy of Cg electrons into account such that the ferromag- 
netic superexchange coupling, which originates from the 
virtue process in Fig. 1(e), occurs, meanwhile the double 
exchange ferromagnetism is compressed completely.Efl 
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FIG. 2. The phase diagram at a; = 0. 



B. Phase separation and A-type antiferromagnetism 

at a; — > 

It was reported experimentally that the phase separa- 
tion was observed in the single crystal of Lai-^Ca^^MnOa 
at a; = 0.05 and 0.08.1 In the present theory, FM coupling 
is very strong near a; = as shown Fig. 2 . Under the 
ferromagnetic background, the AF couplings in He and 
Hyii? are suppressed. The only interacting term survived 
for polarized particles is He -^ in Eq. (|^). Considering 
that the term vanishes unless the neighboring sites are 
occupied by two particles on different orbits at the same 
time, a pure attraction arises between the charge carriers. 
When the lattice distortion increases, FM phase evolves 
into an A-type AF as the distortion forces the orbital bo- 
son to polarize along the c-axis, which further enhances 
the AF coupling via the process (b) . In the case of A-type 
AF, it is FM within the x-y plane, and these FM planes 
are coupled antiferromagnetically. As there is no hop- 
ping between two layers with opposite spins, the system 
can be regarded as a reduced 2-dimensional one. Phase 
separation was discussed in the oiiiij-bajid Kondo lattice 
model numerically and analyticallyll^o and is associated 
with the AF structure. To explore the physical origin of 
PS near a; = 0, the orbital degeneracy of the Cg elec- 
trons has to be taken into account as the PS arises under 
the FM background. The FM near a; ^ in x-y plane 
originates from He ' in Eq. (0), i.e., the spin FM su- 
perexchange process. When the system deviates from the 
point slightly the magnetic structure should not change 
qualitatively. In fact, thfi PS was indeed observed in an 
A-type AF background.Q When the FM forms in the x-y 
plane, the interactions for charge carriers related to AF 
coupling are suppressed, and the attractive interaction in 
the He ' becomes predominant. On the other hand, the 
orbital bosons tend to form an orbital "AF" state, which 
further suppresses the effective hopping term. This prop- 
erty will enhance the relative ratio of the attraction to 
the effective hopping. The strong attraction is the phys- 



ical origin to the phase separation: the charge carriers 
will evolve into two regimes with high and low density 
of the carriers. In the case, the regime with the high 
density of charge carriers has a FM background, and the 
regime with the low density has an AF background as 
only the Jaf term survives at a; = 0. Along the c-axis, 
the antiferromagnetic structure will suppress He and 
He . Thus the pure interaction along the c-axis is repul- 
sive. In reality, the angle between the spins on different 
layer is not absolutely tt. So the dimensionality of A- 
type AF should be between 2 and 3. It is, known that 
a higher dimensionality favors to form PS.E3 The canted 
FM along the c-axis would enhance the stability of PS. 
The phase diagram of PS with respect to the JT effect 
and lattice distortion is shown in Fig. 3. When the lat- 
tice distortion increases, a stronger JT effect is required 
to induce PS_as expected from our general discussion. 
Yunoki et alrB recently reported that the PS appears in 
the phase diagram of a two-dimensional model with the 
Monte-Carlo method. Our result confirms their numeri- 
cal prediction. However, the on-site Coulomb interaction 
was not included in their calculation. 
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FIG. 3. Phase separation and the Jahn- Teller effect near 
a; = 0. 



C. Charge ordering and antiferromagnetism at 

a; = 0.5 

Occurrence of the charge ordering at a; = 0.5 is associ- 
ated with AF structure. For instance, the charge ordering 
with the (tt, tt, 0) pattern occurs under the background of 
C-type AE. that is, AF in the x-y plane and FM along 
the c-axisB The main feature of the charge ordering with 
(tt, tt, 0) is that Mn^+ and Mn'^+ ions aligns regularly in 
the x-y plane and the Mn'^^ ions align along the c-axis. 
It is well known that the CO state is expected to be 
stabilized when the repulsive interaction between charge 
carriers dominates the kinetic energy. The particle-hole 
interaction in the process (b) [or He in Eq. (Q) ] is 
approximately proportional to a;(l — x) and reaches a 
minimal value at x = 0.5. As the sign of the particle- 



hole interaction is negative, it is equivalent to a repul- 
sive interaction between charge carriers (or holes). We 
put forward that the physical origin of the charge or- 
dering results from this process. Of course the direct 
nearest neighbor Coulomb interaction is also favorable 
for the charge. carriers to form CO state, as some authors 
discussed.EilO The phase diagram for charge ordering at 
X = 0.5 shown in Fig. 4 depends on the JT effect and 
lattice distortion as well as the parameters for electronic 
interactions. In the case of C-type AF, the particle- hole 
interaction in He becomes stronger, and the hopping 
term is also suppressed as the Cg electron cannot hop 
to a site with antiparallel spin. Relatively, the effective 
interaction becomes divergent when the spin-spin corre- 
lation becomes AF. In this case the state with a uniform 
density is unstable against the CO. To minimize the po- 
tential energy, the charge carriers tend to form the CO 
within the x-y plane. Along the c-axis, the FM structure 
makes the effective interaction attractive and all charge 
carriers will accumulate along the axis. Therefore, the 
CO has the (tt, tt, 0) pattern. When the JT effect be- 
comes weaker and the lattice distortion increases, the AF 
coupling increases such that the FM coupling along the 
C-axis is suppressed. In the case the effective interaction 
along the c-axis also becomes repulsive. Therefore, the 
G-type AF with the (tt, tt, tt) pattern should be stable. A 
stronger JT effect enhances FM coupling while a stronger 
lattice distortion increases AF coupling along the c-axis. 
Thus it will force the C-type AF to evolve to the A-type 
AF. The charge carriers have the (0, 0, tt) pattern as the 
interaction is repulsive for AF coupling and attractive for 
FM coupling. 
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FIG. 4. The stable magnetic structure at a; = 0.5. The 
G-type AF has a charge carriers distribution with the pattern 
(tt, tt, tt). The C-type AF has the pattern (tt, tt, 0). And the 
A-type AF has a layered pattern (0, 0, tt) 



V. DISCUSSION AND SUMMARY 

Here we discuss the relation between the on-site U and 
sign of effective interactions. The phase diagram in this 



paper bases on the parameters of the model which we 
Ust at the end of Section II. The parameters are roughly 
estimated from the excitation energies of Mn ions and 
the functional density calculations. However the model 
will contain richer phases if we adjust the model param- 
eters. Let us first see two limits, (a) U, Uf ^ JhS, J. 
From the excitation energies of intermediate states we 
list in Table I, the energies in the states of Fig. l(c-d) 
will be much higher than that of the state in Fig. 1(b). 
The FM superexchange coupling will be suppressed, and 
only the AF particle-hole superexchange coupling sur- 
vives. In this case, it favors forming CO in the vicinity 
oi X = 0.5, but does not favor driving charge carriers 
to phase separate near a; = as , equivalently, the net 
particle-particle interaction is repulsive, (b) U ^ JhS. 
The intermediate state in Fig. 1(b) will have a higher 
energy, and the process will be suppressed. The main 
competition comes from the process (c) and (d). The 
former is AF, and the later is FM. Both of them induce 
an attractive interaction for charge carriers. In this case, 
it favors forming PS near x = 0, but does not favor form- 
ing CO at a; = 0.5. The density of dopants of divalent 
ions X in Ri_a;Xa;Mn03 will induce structural parame- 
ters of the crystal, but should have little impact on the 
model parameters related interactions in the ions. Hence 
a mediate value of on-site interaction is very important 
to explain PS and CO simultaneously in the same model 
with the same model parameters. The parameters we 
used in this paper is in the mid-way of the two limits. 
The ratio of U to J^/S is roughly estimated to be 4.8, 
which indicates that the on-site interaction U will also 
have an important impact in the formation of the rich 
phases in the doped manganites. Detailed discussion on 
the model paeameters for transition-metal oxides is re- 
ferred to Ref-CJ 

Next, for doped manganites, the on-site interaction is 
so strong that the system is one of typical strongly cor- 
related electron systems. In our second-order projection 
perturbation approach, part of electron correlations has 
been taken into account, and the direct on-site Coulomb 
interaction and Hund's rule coupling are removed. The 
physical meanings in each term of the effective Hamilto- 
nian in Eq. (0) become much clearer than that in Eq. 
(^. We can see clearly the physical origins of various 
type magnetic structures and related physical processes. 
Strictly speaking, the effective Hamiltonian in Eq. (|^) 
and the original Hamiltonian in Eq. (||) are equivalent 
ONLY in the limit of large U and JhS. Up to the order 
of t/U and t/ JhS, we expect that the two Hamiltonians 
describe the same physics at low temperatures based on 
the principle of the perturbation theory. On the other 
hand, although the projection technique is proved to be 
one of powerfuLtools to deal with strongly correlated elec- 
tron systems,c3 it is still a challenging problem to deal 
with the effective Hamiltonian. Our spin-orbital-charge 
coherent state theory is just an initial step to understand 
the physics in doped manganites. 

Finally, we come to discuss the relation of the charge 



inhomogeneity and the long-range Coulomb interaction. 
The tendency to PS and CO depends on the magnetic 
structures. It is anticipated that the inclusion of the 
longer-range Coulomb interaction will lead— to a sta- 
ble and microscopically inhomogeneous state.L3 Recently 
Mori et al. reported that charge stripes arise iit-jthe 
range of a; > 0.5, and tend to form stripe pairsJlJ It 
reveals strong repulsion for the charge carriers since the 
striped Mn'^+ ions are separated by Mn*+. By consid- 
ering the nearest neighbor interaction and longer-range 
Coulomb interaction. Low et al. obtained some numer- 
ical evidences to support ptiie stripe phase at x=0.5 in 
a simplified lattice modelMs Our electronic model does 
not include the long-range Coulomb interaction. How- 
ever even if we take the interaction into account, the 
nearest neighbor interaction should not be very strong 
such that it does not destroy the tendency to PS in the 
vicinity of a; = caused by the superexchange attraction. 
Of course, it enhances tbe [tendency to CO near x — 0.5 
as some authors argued. O'Eil 

In short, starting from a realistic model, we derive an 
effective Hamiltonian by means of the second order pro- 
jection perturbation approach in the case of the strong 
Hund coupling. In order to treat the model, we intro- 
duce a new type of bosons for orbital degree of freedom 
as well as bosons for spins and fermions for charge car- 
riers. A spin-charge-orbital coherent state theory is pre- 
sented. Physically, by adjusting the orbital ordering of 
the charge carriers, we find that the JT effect and the lat- 
tice distortion have strong impact on the electronic col- 
lective behaviors as well as the magnetic structures. At 
the undoped case (a; = 0), the ferromagnetism originates 
from the FM superexchange coupling, and the anisotropy 
of A-type AF is induced by the crystal field. Away from 
the undoped case, the FM superexchange coupling term 
is responsible to form PS near the slight doping regime 
while the AF particle-hole interaction drives the charge 
carriers to form CO near x = 0.5. 

This work was supported by a CRCG research grant 
at the University of Hong Kong. 



APPENDIX A: THE SECOND-ORDER 
PERTURBATION EXPANSION 

In this Appendix, we derive the effective Hamiltonian 
in Eq.(0) in which correction of the finite but large Jh 
effect is taken iato account. We follow Schrieffer and 
Wolff's methocCa to derive the Kondo HamiltonijaB from 
the periodic Anderson Hamiltonian (Also see Ref .c2l) . Ac- 
cording to the projection operator P, the Hilbcrt space 
is divided into a subspace P, which consists of holes and 
single occupancies with spin 5* + 1/2, and a subspace 
Q{= 1 — P) with at least one double occupancy or one 
single occupancy with S— 1/2. The Schrodinger equation 
for the system is written as 



H\ 



= E„ 



where Eg is the ground state energy. The equation can 
be expressed in the two subspaces P and Q 



PHP\4>)+PHQ\4>) ^EgP\ 
QHP\4>)+QHQ\4)) ^EgQ\ 



(Al) 
(A2) 



The Haniihonians PHP and QHQ act within the sub- 
spaces P and Q, respectively. PHQ and QHP con- 
nect the two subspaces. To ehminate the state Q |<^) in 
Eq.(Al), we reduce the problem in the subspace P. 



[H, 



E„)P\<h) =0 



where 



Hp = PHP - PHQ- 



1 



QHQ - Eg 
The operator Q can be expanded as 



QHP. 



In the large U one-band Hubbard model, we have only 
one projection operator for the intermediate state, i.e., 
the operator for the double occupancy. In our case, we 
consider four intermediate states. The four energy differ- 
ences and projection operator Qi^a of the intermediate 
state shown in Fig.l (b-e) are listed in the Table I. 

TABLE I. The four energy differences and projection op- 
erator Qi,a of the intermediate state shown in Fig. 1 (b-e). 



Energy AEa 



Jh{2S + 1) 



U + JhS 



U' + P+.JhS 



U'-hJ 



Operator Qi, 



, r' r- , P 

2S + 1 /n-n-' '.7,°- '■T-°' '■' 



^i,'y,| '''ij'YiJ, 



-Si-a + (S + l/2)l \ ^t 



25+2 ; ''i,7,CT'^J,7,o-'-'».7.3 
njzl 



Si-a + (S + 3/2)l \ „t „ ,p 

2S+2 j S,7,ct''1,7.o- -■ ^■'<:'> 



Eq^,AyI(p^,+ Y.p^ 



+ ) 



Ji^i 



7=±1 






where Qa is the projection operator that there is at 
least one double occupancy or single occupancy with spin 
25—1. For our purpose, we just consider the energy cor- 
rection up to the second order perturbation. Hence we 
take approximately 

i j¥'i 7=±1 

One of important properties is 

As we are merely concerned with the low-energy exci- 
tation, the term 



^ a ' 



1 



QHQ - Eg 
is replaced approximately by 



Qo 



AEc 



where AEa is the energy difference of energy with one Qa 
and energy of PHP. Thus the Hamiltonian is reduced 
to 



H. 



eff 



PHP - J2 -^PHQaHP 



A^a 



z. 

K. 



1 

2 ■ 

5 A. 
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